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L  Rttiarcli  AcoqnipliilnaMis. 

FM  vorla  iMllwdt. 

We  developed  a  fast  vonex  method  for  solving  the  incompressible  Euler  etpntktts  in 
three  dfanoisions.  It  is  based  on  a  combination  of  Andenon's  Method  of  Loosl  Gonee- 
tkms,  wMdi  uses  a  hybrid  of  a  finite  tfifferenoe  repiesentation  and  a  paniele  repiesen- 
tadon  10  oonqmte  the  velocity  field  induced  on  foe  vortices,  and  of  adaptive  mesh 
lefinement  for  the  finite  difference  calculation.  The  resulting  algorithm  is  foster  than 
any  other  existing  finite  difference  or  mulfipole  •  based  algori^  for  this  problem,  and 
is  10  -  20  times  foster  than  the  direct  N*body  method  for  problems  of  current  research 
interest.  In  addition,  we  developed  a  good  deal  of  technology  for  finiie*difforence 
Poissmi  stivers  of  independent  interest  and  usefulness,  such  as  a  muldgrid-based 
ktcal-iefiaement  algorithm  for  Mehrstellen  discretizations  of  Poisson's  equation  in 
three  dimensions,  and  efficient  treatment  of  boundary  conditions  for  boundaries  that 
have  infinite  extent  in  one  or  more  of  foe  coordinate  dfaeetions. 

Projeclioo  mellKNbi 

We  continued  our  work  on  projection  methods  for  the  inconqxtssible  Euler  and 
Navier-Stokes  equations  in  two  and  three  dimensions.  These  m^x)ds  use  the  hi^ 
resolution  upwind  finite  difference  methods  developed  for  hypeibdic  problems  to 
evaluate  the  nonlinear  advective  terms,  combined  with  a  second  order  predictor* 
corrector  time  differencing  to  intertwine  the  projection  operator  and  the  viscous  terms. 
We  have  developed  new  veiiions  of  the  second*order  algorithm  of  Bell,  Colella,  and 
Glaz  that  are  more  dlicient  than  the  previous  versions,  increasing  the  largest  peimissi> 
bk  time  step  by  a  factor  of  two,  and  systematically  introducing  multigrid  tinoughout 
for  solving  the  linear  systems  arising  tmn  the  discretizations  of  the  elliptic  and  para* 
bolk  (q^iors.  Overall,  the  algorithm  is  a  factor  of  five  faster  than  the  previous 
method.  We  have  extended  titis  version  of  the  algorithm  to  ^neral  variable-density 
incompiessibk  flows  and  to  three  tfimensions.  We  have  also  extended  the  projection 
formalism  to  9pp\y  to  the  zeto-Mach  number  model  for  reacting  fluid  flow  propt^  by 
Majda  and  Sethian.  Using  this  approach,  we  are  able  to  represent  the  effects  of  large 
volumetric  expanrions  of  the  fluid  due  to  the  energy  release,  and  voiticity  produc¬ 
tion  due  to  large  density  variations,  without  an  explicit  treatment  of  acous^  waves. 
This  leads  to  an  efficient  computational  method,  since  our  approach  can  use  a  time 
step  that  is  limited  by  the  CFL  condition  for  the  advective  transpoit,  rather  than  that  of 
the  sound  waves,  as  would  be  the  case  in  a  fully  compressible  calculation. 

Finit#  diffortnce  methods  for  hyperbolic  probkitts. 

Wc  have  continued  our  woifc  in  high-resolution  finite  difference  methods  for  hyper¬ 
bolic  equations  based  on  second  onler  accurate  extensions  to  Godunov's  method.  We 
have  d^loped  an  iiiq>licit/bxplicit  method  for  hypeibolk  and  hyperbolic/lpanbolic 
proMems  in  one  apace  dimension,  that  has  the  property  that  the  finite  difference  stencil 
switch  ftom  being  explidl  to  implicit  in  time  continuously  and  locally  in  each  cell  and 
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for  each  chancteritic  family,  (fepending  on  whether  the  local  CFL  number  for  that 
characteristic  is  smaller  or  larger  than  one.  In  the  limit  where  all  of  the  CFL  numbers 
are  smaller  than  one,  the  iiie£od  is  a  standard  second  order  finite  difference  method; 
to  the  limit  of  all  ehaiicteristics  in^>lieit  and  the  solution  at  steady'State,  the  method  is 
a  seeood«oider  accurate  discretna^  for  the  neady  state  equadons.  We  have  also 
extended  the  second-order  Oodimov  method  of  Bell,  Colella,  and  Ttongenstein  for  non- 
convex,  non-striedy  hyperbolie  systems  to  the  equations  of  ideal  MHD  in  one,  two  and 
three  dimensions.  This  system  exhibits  a  much  richer  variety  of  wave  propagation 
behavior,  which  had  to  be  analysed  in  order  for  these  methods  to  be  extended  to  that 
case,  to  more  than  one  dimension,  we  use  a  {uqjectiott  operator  similar  to  that 
developed  for  toconqxessible  flow  to  insure  (hat  the  divergence-ftoe  constraint  mi  the 
magnetic  field  is  satisfied.  Finally,  we  developed  a  mudi-simplified  verdon  of  the 
secondmrder  Godunov  method  for  general  equations  of  state,  in  the  context  of  the 
multUhiid  algorithm  discussed  below.  As  well  as  being  a  factor  of  three  faster  than  the 
previous  algorithm  of  Oolella  and  Glaz,  it  is  much  more  general,  qiffiying  to  the  case 
of  equations  of  state  for  which  neither  the  pressure  nor  internal  energy  necessarily 
remain  positive,  such  as  arise  to  representing  liquids  and  compressive  modons  to 

Adaptive  nclhoda. 

Our  primary  emphasis  has  been  to  the  development  of  tracking  methods  for  represent¬ 
ing  seketed  sharp  fronts,  toterfhces,  and  complex  boundary  geometries  for  conquessl- 
ble  flow  problems.  Our  approach  is  based  on  uring  a  volume^-fiuid  descnpiim  for 
the  tracked  fiont,  while  capturing  the  remaining  flow  features  using  a  conservadve 
finite  (Ufference  sriieme  on  either  skle  of  tracked  firont.  For  the  case  of  tracking 
ihoeks  and  ctmiidex  solid  bodies,  we  deal  with  the  stability  problems  arising  fiom  con¬ 
servative  finite  (fifferencing  on  small  cell  fragments  by  uring  the  flux  icdistiibudon 
scheme  of  Oiem  and  Colella.  We  develqied  such  an  algorithm  for  the  case  of  an 
unsteady  tracked  shock  to  two  dimensions,  for  whidi  the  volumcKiMluid  representa¬ 
tion  allowed  us  to  conecdy  represent  a  shock  tindeigoing  larp  disiordons  and  changes 
to  topology.  We  have  coupled  that  method  to  the  block-structured  adaptive  mesh 
refinement  algorithm  of  Berger  and  Golelia  to  a  comideiBly  general  fashion  that  allows 
the  tracked  fiont  to  cross  the  boundaries  between  different  refinement  levels.  We  have 
also  used  the  same  code  structure  to  develop  a  Cartesian  mesh  representation  of  oom- 
plex  boundariea,  to  which  the  boundary  the  body  is  treated  using  the  same  tech- 

idques  as  the  ttaeked  from.  The  modification  of  the  two  dimeniioaal  trackiag  code  to _ 

do  the  Canerian  grid  case  was  routine,  taking  about  two  weeks;  the  threeHUmensional  For  y 
veffioo  of  the  Gvterian  grid  geometry  algorithm  is  in  the  testing  sta^  Rnally,  we  a&i 
have  developed  a  version  of  ite  nxire  classical  muhifluld  algorithms  for  representing  ^  □ 

the  Inteifhce  between  two  or  more  diflerent  materials.  The  prima^  innovation  to  this  ^ 

work  is  toil  our  algorithm  is  dieimodynamicaUy  consistent,  treating  materials  to  coo-  • 

tael  with  one  another  having  widely  varying  cmxq^tessilnUtiet  in  a  fashion  consisiem 
with  the  hypoiheris  that  the  pressure  to  conttouous  across  the  interface.  Our  algorithm 

can  be  formaDy  derived  as  a  discrefizaUon  a  system  of  partial  diffbtential  equatioos  ^ - 

for  a  tiioltioom|N»eot  titixture  that  ean  be  proven  to  be  hypeihttiic  to  the  case  of  — 

smooth  promne  variations,  with  the  condition  of  pressure  equilibrium  among  the  Special 
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componemt  at  eveiy  point  in  q>aee  hiding  as  an  initial  value  eonitraint  !n  cases 
where  large  an^titude  pressure  waves  cause  the  oontinuity  of  the  pressure  hypothesis 
to  be  vkdaied,  we  introduce  a  relaxation  scheme  to  restore  {vessure  equilibrium  among 
the  fluid  components.  We  have  also  coupled  this  algorithm  to  the  adigxive  mesh 
methods,  which  hu  proved  panicttlariy  useful  in  our  investigations  of  shocks  interact* 
ing  with  interfroes  between  two  gases. 

Appiicatioiia. 

We  applied  the  numerical  methods  described  ri)ove  in  a  variety  of  spocific  physical 
apfdii^ons,  mainly  in  the  area  of  shock  dynunics.  We  have  continued  our  work 
planar  shock  n^ec^  and  reftaction.  We  used  a  shock  tracking  /  adtqnive  refinement 
code,  along  with  carefully  derigned  shock  tube  experiments,  to  determine  the  criterion 
for  tranrition  between  regular  and  Mach  reflection  for  strong  shocks,  in  a  reghne 
where  the  analytical  theory  allows  both  solutkms.  We  found  that,  fbr  inviscid  prob¬ 
lems,  Mach  reflection  is  the  correct  soludon  in  die  indeterminate  region  of  parameter 
qpace,  wifli  viscosity  cansing  a  displacement  of  the  transition  point  towards  earlier  reg¬ 
ular  reflection.  In  our  work  on  refiicdon,  we  investigated  the  collision  of  a  planar 
shock  with  an  oblique  planar  interface,  where  the  speed  of  sound  on  the  side  of  the 
interface  qpposite  the  shock  was  lower  than  that  of  the  initially  shocked  gas.  The 
numerical  method  used  was  the  multifluid  method  described  above,  coupled  to  the 
atUqttive  mesh  refinement  scheme.  Bnally,  we  investigated  the  use  of  our  high- 
reatdudon  finite  difference  methodology  for  the  case  of  flow  in  parous  media  widi  real¬ 
istic  heterogeneities.  We  used  $tate-of-the-an  geostadsdcal  techniques  to  generate  per¬ 
meability  and  porosity  fields  on  rectangular  grids  from  well  data  from  the  Wilmingion 
oil  field  in  Long  Beach.  CA  .  We  found  that  apjarquiately  derigned  versions  of  our 
finite  difference  methods  for  the  transport  equatUms,  plus  a  suitable  variation  on  mol- 
dgrid  for  the  pressure  equation,  leads  to  an  efficient  and  accurate  treatinem  of  the  flow 
in  such  fields. 
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